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We address the question: How does capillarity propel microspheres along curvature gradients? 
For a particle on a fluid interface, there are two conditions that can apply at the three phase 
contact line: Either the contact line adopts an equilibrium contact angle, or it can be pinned by 
kinetic trapping, e.g. at chemical heterogeneities, asperities or other pinning sites on the particle 
surface. We formulate the curvature capillary energy for both scenarios for particles smaller than 
the capillary length and far from any pinning boundaries. The scale and range of the distortion 
made by the particle are set by the particle radius; we use singular perturbation methods to find the 
distortions and to rigorously evaluate the associated capillary energies. For particles with equilibrium 
contact angles, contrary to the literature, we find that the capillary energy is negligible, with the 
first contribution bounded to fourth order in the product of the particle radius and the deviatoric 
curvature. For pinned contact lines, we find curvature capillary energies that are finite, with a 
functional form investigated previously by us for disks and microcylinders on curved interfaces. 
In experiments, we show microsphere migrate along deterministic trajectories toward regions of 
maximum deviatoric curvature with curvature capillary energies ranging from 6x10^ — 5x10^ ksT. 
These data agree with the curvature capillary energy for the case of pinned contact lines. The 
underlying physics of this migration is a coupling of the interface deviatoric curvature with the 
quadrupolar mode of nanometric disturbances in the interface owing to the particle’s contact line 
undulations. This work is an example of the major implications of nanometric roughness and contact 
line pinning for colloidal dynamics. 


I. INTRODUCTION 

Capillary interactions are ubiquitous between particles 
on fluid interfaces. They trap particles at the interface 
[T] and determine their ensuing organization l2]-[4], 
allowing particles to be widely exploited in techno¬ 
logical applications such as stabilization of foams and 
emulsions [5], and in settings as diverse as the food [ 6 ], 
pharmaceutical [7], mineral recovery [ 8 ], and petroleum 
industries [9]. For microparticles of radius a at planar 
interfaces of tension 7 , gravity is irrelevant, as the Bond 
number Bo = <C 1, where g is the gravitational 

acceleration constant, and Ap is the density difference 
between the two fluids on either side of the interface. In 
this limit, particles with undulated contact lines distort 
the interface around them, with an associated excess 
interfacial area. The deformation fields and interfacial 
area depend on the relative position of the particles, 
yielding decreasing capillary energy as particles approach 

m- 

In this research we are interested in the behavior of 
isolated microparticles trapped on curved interfaces. 
In experiment, microparticles migrate along curvature 
gradients to sites of high curvature, as has now been 
observed for microcylinders [ 2 ], microspheres m, and 
microdisks m- Theoretically, the curvature capillary 
energy driving this migration is simply the sum of the 
surface energies and pressure work for particles at the 
interface. When particles attach to their host interfaces, 
they change the interface shape owing to the boundary 
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condition at the contact line where the interface meets 
the particle. There are two limits for this boundary 
condition. The interface can intersect the particle 
with an equilibrium contact angle Oq, determined by 
the balance of surface energies according to Young’s 
equation m- Alternatively, the contact line can be 
pinned by kinetic trapping at heterogeneities, roughness 
or other pinning sites on the particle surface [iniiiiHiH]. 
Curvature capillary migration depends on the coupling 
of the particle-sourced distortion with the host interface 
shape. For a particle which is much smaller than the 
capillary length, far from pinning boundaries, and there¬ 
fore small compared to the prevailing curvature fields, 
particle-sourced distortions decay over distances compa¬ 
rable to the particle radius. This scenario lends itself 
to analysis by a singular perturbation method [ 20 ] . 
In a particle-fixed reference frame defined in the plane 
tangent to the interface, the host interface shape 
around the particle can be locally decomposed into 
terms proportional to the host interface mean curvature 
and deviatoric curvature. By rescaling the governing 
equations with respect to the particle radius, the shape 
of the interface around the particle can be found over 
length scales comparable to this radius. This approach 
allows clear treatment of limiting boundary conditions 
on particle sourced distortions and the ranges of validity 
of local expansions of the host interfaces. It also clarifies 
the appropriate outer limits on (area or contour) inte¬ 
grals used to evaluate the excess area of the interface 
created by the particle. Using this method, we have 
previously solved the curvature capillary energy for 
disks with pinned contact lines, which applies equally to 
spheres with pinned contact lines in the limit of weak 
undulations [ 12 ]. We apply the approach to derive the 
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FIG. 1. Schematic representation of a sphere trapped at a 
curved interface and its mapping to the tangent plane. 


curvature capillary energy for spheres with equilibrium 
contact angles. While the curvature capillary energies 
for this scenario have been derived previously and 
reported to be quadratic in the deviatoric curvature 
of the interface [niiiiHis], we find that this term has 
prefactor zero. We identify the source of the discrepancy 
between our result and that published previously. We 
predict that spheres migrate by capillarity on constant 
mean curvature interfaces only if their contact lines are 
pinned, and that spheres with equilibrium contact lines 
would have energies several orders of magnitude weaker 
than is observed in experiment. 

We perform experiments in which we record the 
trajectories of polystyrene microspheres at hexadecane- 
water interfaces with well defined curvature fields. We 
compare the energy dissipated by particle migration to 
the curvature capillary energy expressions, and find that 
the spheres migrate in agreement with the expression 
for pinned contact lines. 


II. THEORY 


Here we give a concise derivation of the capillary en¬ 
ergy of a sphere trapped on a curved fluid interface. 
Without loss of generality, we focus on interfaces with 
zero mean curvature as the role of finite mean curvature 
gradient has been addressed in the literature [T2j [22] . 

We first consider the free energy of a system including 
a free interface with a spherical colloidal particle sub¬ 
merged in subphase fluid 1 . In this case, the free energy 
can be written as. 


El = 7 


+ ( 1 ) 


where is the Monge representation of the host interface 
height prior to particle attachment, D denotes the entire 
interfacial domain, dA = dxdy, 71 is the surface energy 
of the colloid with fluid 1, and As = 47ra^ accounts for 
the total area of the colloid. In writing the above expres¬ 
sion we have assumed small slopes, i.e. |V/io| ^ 1 - Upon 
attaching to the host interface, the particle creates a dis¬ 
turbance to satisfy its boundary condition on the three 
phase contact line. The free energy of the system in this 
case is, 

E 2 = 71^1 + 72^2 + 7^(1 + ^A^)dA, ( 2 ) 

D-P 

where h is the interface height after the particle is ad¬ 
sorbed, 7 iHi and 72^2 are the product of the surface 
energies and wetted areas for the solid in contact with 
the upper and lower fluids, respectively, and P denotes 
the domain under the particle. The capillary energy as¬ 
sociated with adsorption of a particle to an interface can 
be found simply via subtraction of the two energies, 

E = E2- Ei= 7iHi + 72H2 - 71 + 

D-P P 

( 3 ) 

where 7 is the disturbance field defined as 7 = h — ho. 
In the above expression, the first integral is due to the 
disturbance created by the particle and the second inte¬ 
gral is the area of the hole in the host interface under 
the particle. We integrate the first integral by parts and 
then apply the divergence theorem to have, 

^ [ Vr; Wt? ^ ^ ^ho]dA = 

D-P 

j ) + 7V/10] • 

d{D-p) 

( 4 ) 

where d{D — P) denotes the contours enclosing the do¬ 
main D — P (See Fig. 1). There are two contours enclos¬ 
ing this domain: one, not shown, infinitely far from the 
particle, with outward pointing vector in the radial direc¬ 
tion, and the other enclosing the region P with outward¬ 
pointing unit normal vector m. Consequently the cur¬ 
vature capillary energy of the system can be expressed 
as, 

E = E2 — El = jiAi + 72^2 — + 

7 [ ^ [^ri\7ri + riVho]-mds-^{l + ^EEE^dA]. 

d(D-P) P 

( 5 ) 

We expand the interface locally around an arbitrary 
point to obtain a saddle shape, 

ho{x,y) = -y^), 


D 


2 


(6) 
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FIG. 2. Schematic representation of a sphere trapped at an 
interface with (a) a pinned undulated contact line and (b) an 
equilibrium wetting boundary condition with a contact angel 
Oo. 


where (x, y) coordinate is tangent to the host interface, 
oriented along the principle curvatures ci and C2, and Acq 
is the deviatoric curvature of the host interface defined 
as, 

Aco = Cl - C2, ( 7 ) 

where we adopt a convention such that ci is always posi¬ 
tive. We discuss the range of validity of Eq. |^in terms of 
an asymptotic treatment in Appendix A, and in terms of 
a comparison of host interface shape used in our exper¬ 
iments in Fig. To evaluate the integrals in Eq. we 
must define the boundary condition at the three phase 
contact line, and determine the associated disturbance 
field T]. We discuss two distinct scenarios for this energy 
(see Fig. [^. 


A. Pinned contact line 

The curvature capillary energy E, for a particle with 
nearly circular cross section and a pinned contact line 
trapped on a host interface with arbitrary mean curva¬ 
ture Hq and deviatoric curvature Acq was derived pre¬ 
viously m- The height of the pinned contact line con¬ 
tour can be decomposed into a multipole expansion with 
quadrupolar mode of amplitude hp. The associated cur¬ 
vature capillary energy is. 



FIG. 3. (a) Schematic view of a sphere of radius a at a liquid 
interface in x — z plane, where Oq and ro are the contact 
angle and the radius of the contact line for a flat interface. 
ho describes a curved interface without a particle, h = ho -\-r] 
accounts for the height of interface in presence of the particle. 
0 is the polar angle of the deformed contact line, and p is its 
radial position, (b) Schematic view of the contact lines in 
X — y plane in flat (in blue) and curved interfaces (in red). 


B. Equilibrium wetting 

In Appendix A, we use a formal matched asymptotic 
expansion [niEn] to evaluate the curvature capillary en¬ 
ergy in terms of the small parameter A = roAco, where 
ro = a sin6>o is the radius of the hole made a sphere is 
a planar interface. Here, we give the main features of 
this derivation in dimensional form, and note the bounds 
on the solutions and ranges of validity implied by the 
asymptotics. 

For a spherical particle in a curved interface with an 
equilibrium wetting condition, the fluid interface deforms 
until it satisfies the equilibrium contact angle at every 
point on the contact line. The contact line shape is not 
known a priori and must be determined as a part of anal¬ 
ysis, as was originally done by Wiirger m- The shape 
of contact line can be deduced from geometrical relation¬ 
ships (see Fig. to be. 


^ ^ h(r e p) 

cos U = cos dp H-, 

( 9 ) 

a 

sin 0 = -, 
a 

(10) 


E = Eo- 


77ra^( 


hpAcp 

2 


4 


( 8 ) 


In this expression, the first term Ep is independent of the 
local curvature. The second term predicts that a parti¬ 
cle will move to sites of high deviatoric curvature, while 
the third predicts particle migration along gradients of 
mean curvature. To understand the relevant importance 

of these terms, we consider , ^ ° ^ where Ap is 

the pressure jump across the interface (see Appendix C 
for non-dimensionalization scheme). This ratio suggests 
that, for sufficiently small pressure jump across the in¬ 
terface, the effect of mean curvature can be neglected. In 
section 3 , we explore this regime in experiment. 


where 6> is a polar angle which located the contact line in 
a spherical coordinate system with respect to the particle 
center, cos^o = ^5 and p{(j)) is the contour of the contact 
line projected into a plane, which obeys, 

p = a \/1 — cos^^. (11) 

By substituting Eq.[^and[TQ|in this expression and as¬ 
suming ^ 0 (A) <C 1, the radial location of contact 

line projected into the x — y plane is, 

1 

p = ro- ft.(r)cot6»o - -- -{cot'^do + 1) + --Irep) 

z rp 

( 12 ) 
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The corresponding unit vector m and the arc length 
element ds up to the first correction in A are, 


m = — — 


cot Oq dh{r) 
To dcj) 


IrGp’ 


ds = ro[l - cot6>o^^]#|^^^. 


^0 


(13) 

(14) 


Young’s equation requires, 


np • n/ = cos6>o|rep. 

where np is the unit normal to the particle and n/ is 
the unit normal to the interface evaluated at the contact 
line. By applying this expression (see Appendix A for 
details), the boundary condition for the interface shape 
at the three phase contact line to the leading order is. 


h dh 
— = ^lr=ro- 


(16) 


The height of interface satisfies the Young-Laplace 
equation, which, assuming small slopes, reduces to the 
Laplacian: 


= 0 . 


(17) 


agrees well with the full function describing the inter¬ 
face shape /i,micropost. quadrupolar mode of /^^rnicropost 

agrees with Eq. to better than 2% error. This is the 
sole mode that can couple to the quadrupolar undulation 
of the particle. While we perform this matching in de¬ 
tail only for the host interface used in our experiments, 
the asymptotic analysis can be applied similarly to any 
fluid interface provided ^ ^ 1, Rc is the characteris¬ 
tic length of the host interface. The deviatoric curvature 
part of far field boundary condition is given by Eq. 
with the details for the curvature source (e.g. the micro¬ 
post in our example) emerging in the expression for Acq. 
In the event that the interface has finite mean curvature, 
an additional term emerges with the local mean curva¬ 
ture of the interface. 

We evaluate 

(i) the self-energy of the disturbance created in the host 
interface, the sum of the following two contour integrals: 



p 



• mds =ro 



Acocos^20 

144 


d(p 


144 


( 21 ) 



• mds =0 


r^oo 


( 22 ) 


The boundary condition for the particle-sourced dis¬ 
tortion in the far field can be derived using a Van Dyke 
matching scheme [H [20] to match the interface shape 
far from the particle to the host interface shape. The 
resulting boundary condition far from the particle is, 

lim h(r, 0) = ho = cos 20, (18) 

r^oo 4 


where the limit r ^ oc implies exploring regions around 
the particle large compared to the particle radius. Using 
Eq. \T6\ and Eq. [^the leading order interface shape is, 

, , Acocos20 , 9 rn . 

Hr,^)= % (19) 

The corresponding disturbance to the interface is. 


T] = h — ho = 


Aco cos 20 Vq 
12 


( 20 ) 


This term captures the increase in interfacial area owing 
to the particle sourced distortion. 

(ii) the interaction of the disturbance and the host in¬ 
terface, given by the sum of the following two contour 
integrals: 



24 

(23) 



r^oo 


ijS/ho • mds =ro 



Acocos^20 

24 


dcj) — + 


24 


(24) 


These terms are equal and opposite, and hence sum to 
zero. Thus, the net contribution of the particle induced 
disturbance, the sum of (i) and (ii), is 


Asymptotic analysis for the curved interface around the 


circular microcylinder (Appendix A) show that Eq. 18 


mains valid everywhere the disturbance (Eq. 20) is finite. 
With this information, we determine the curvature capil¬ 
lary energy Eq. In so doing, we can straight-forwardly 
apply the limits at the contact line and r ^ oc, know¬ 
ing that these functions remain valid within this domain. 
This is indeed consistent with experiment. We mold the 
fluid interface by pinning it to a micropost. The interface 
forms a shape with a well defined curvature field. Eor 
micro-particles on interfaces with deviatoric curvatures 
like those in our experiments Acq ^ 400 — 3000 
the distortions will be <C 10 “within 5 — 6 radii, a, 
from the the center of particle. Over such ranges, as 
shown in Eig. the shape of the fluid interface Eq. 


/ -I A 2 4 

[-?7V^ + 7?Vfeo] ■mds= (25) 

d(D-P) 

(iii) We calculate the energy decrease owing to the area 
of the hole under the particle as. 





-nrl 


P 


144 


(26) 


The second term in the above expression is equal and 
opposite to the term in Eq. Summing the curvature 
dependent terms in (i),(ii) and (iii), the energy costs of 
increased area by the deformation field is shown to be off¬ 
set by the energy decrease owing to the area eliminated 
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FIG. 4. Comparison of the shape of the host interface used 
in experiment, graphed with respect to the plane at an arbi¬ 
trary point L = 400 (im from the micropost center /^^"^icropost 
(Eq. 1^ (right) and Eq. the limiting form of the inter¬ 
face far from the particle lim /i(r, 0) (left). The quadrupolar 


mode of 


agreees with Eq. 18 to within 2% error. 



under the particle. This exact canceling of two effects 
occurs both for the pinned and equilibrium contact lines. 

To complete this derivation, the contribution of the 
(curvature independent) wetting energies can be simpli¬ 
fied by noting, 


71^41 = 7i27ra^(l + cosOq), 

(27) 

72^2 = 7227ra^(l - cos^o), 

(28) 

7i^s = 7i47ra^, 

(29) 

a 72-71 

COSC/q = -, 

7 

(30) 


summing these contributions, the net curvature capillary 
energy to order is identically zero, i.e.. 


E 

77rrg 


Ej, + 0{X^), 


(31) 


where the first term in this expression is Pieranski’s trap¬ 
ping energy [T] for a sphere obeying an equilibrium con¬ 
tact angle in a planar interface Ep = — sin“^6>o(l — 
cos^o)^- The above expression is exact up to the 
term; the determination of this bound is described in 
Appendix A. 

Eq.[3T] shows that there is no quadratic term in Acq for 
spheres with equilibrium contact lines, and that the cur¬ 
vature capillary energies for these particles are extremely 
weak. A similar conclusion has been found for perfect 
disks with circular pinned contact lines [12]. These re¬ 
sults differ significantly from prior theory in the litera¬ 
ture for this problem. The origin of the discrepancy is an 
inappropriate treatment of the contour integral given in 
Eq. which was assumed to be zero in prior work. In 
Appendix B this discrepancy is explored in greater detail 
by evaluating the integrals in questions three ways, all of 
which agree with our result. We also follow the approach 
in the literature for two forms of the contour integral for¬ 
mulation that are, in principle, equivalent, and show that 
they lead to contradictory results. 


EIG. 5. (a) Schematic view of curvature field around the mi¬ 
cropost. (b) SEM image of a microsphere and PDMS negative 
of the air-water interface showing particle roughness and (G) 
AEM reveals nanoscopic roughness of the sphere surfaces with 
RMS values ranging of 15 — 21 nm. 


III. EXPERIMENTS 


We study migrations of polystyrene colloidal spheres 
(Polyscicences, Inc.) with mean diameter of 2a = 10 fim. 
Fig. [^(b) illustrates an SEM image of a microsphere re¬ 
vealing qualitatively the surface roughness of the particle. 
AFM measurement (Bruker Icon) indicates that the root 
mean squared roughness of the particles is ^ 15 — 21 nm 
(see Fig.|^(c)). 

We impose a curvature field to the host interface using 
a technique reported previously [2| which we recapitulate 
briefly. A curved oil-water interface is formed around a 
micropost which is either circular or square in cross sec¬ 
tion (see the schematic in Fig. |^a)). The interface pins 
to the edge of the post, and has a height Hm at the post’s 
edge. The post is centered within a confining ring located 
several capillary lengths away = 5.5, where Vring 


is the radius of the outer ring. By adjusting the volume 
of water, the slope of the interface at the post’s edge is 
adjusted to be t/; ^ 15 — 18°. This system is gently cov¬ 
ered in hexadecane in order to prevent evaporation and 
to protect the interface from stray convection. The in¬ 
terface height in a region sufficiently close to the circular 
post is well approximated by: 


/^micropost = tan ^ ln( ^), (32) 

where L is the distance to the center of the post. This 
interface has zero mean curvature Hq = (ci + C 2)/2 = 0, 
and finite deviatoric curvature Acq = ci — C 2 varying with 
the radial position. Owing to the finite volume of fluid, 
there is a weak but negligible pressure drop across the in¬ 
terface as we confirm in numerical simulations described 
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FIG. 6. Trajectories of microspheres around a microcylinder, 
(a) Time stamped image of microspheres trajectory (scale bar 
= 100 /xm). (b) Trajectories for spheres migrating on the in¬ 
terface. Colors for trajectories in (b) correspond to trajecto¬ 
ries of similar colors in (a). 


in Appendix C. 

A dilute suspension of microspheres in hexadecane is 
prepared. A drop of this suspension is carefully dropped 
on top of the oil phase. The particles then gently sedi¬ 
ment under gravity. Once attached to the interface, they 
migrate uphill in a deterministic path along deviatoric 
curvature gradients. We only focus on isolated particles 
far from neighbors (distances greater than 10-15 radii) 
and the micropost {L > Rm + 10a) to rule out the pair 
capillary and hydrodynamic interactions M- The capil¬ 
lary energy was estimated by evaluating the total dissipa¬ 
tion according to the appropriate drag formula (Stokes’ 
law) along particle trajectories. 


IV. RESULTS AND DISCUSSIONS 

In Fig. I^a), we illustrate the time stamped images of 
trajectories for migrating spheres for constant time incre¬ 
ment (At = 1 s). These images reveal that the spheres 
are propelled faster in the region closer to the post where 
the magnitude of deviatoric curvature is greater. Note 
that the size of the spheres are so small that the inertial 
effects can be neglected {Re ^ 10“^) within the entire 
trajectory. These trajectories are nearly radial, as the 
corresponding curvature field around the cylindrical post 
has no dependency on the azimuthal angle (j). Fig. [^b) 
shows the radial distance of the migrating microspheres 
from the center of the post, L, as a function of time re¬ 
maining until contact, t —tc, where tc is the time in which 
the sphere reached the edge of the post. Qualitatively, 
these trajectories are similar to those reported previously 
for microdisks with pinned contact lines migrating in cur¬ 
vature fields. To investigate this quantitatively, we com¬ 
pare energy dissipated along the particle trajectory to 
theory. 

To do so, we note that, in the limit of zero inertia, 
and neglecting potential energy differences, the curva¬ 
ture capillary energy expended to drive the particles is 


xIO'® 



FIG. 7. Gomparison of predicted curvature capillary energy 
for particles with pinned contact lines (Eq. (solid lines) 
against those extracted from experiment for trajectories in 
Fig. I^b) for isolated particles (symbols). The colors corre¬ 
spond to the trajectories in Fig. HJb). The linear fit is ex¬ 
cellent for all trajectories, with coefficient of linear regression 
R^ — 0.999 for the worst case. 


balanced by viscous dissipation. The total energy dissi¬ 
pated along the trajectories can be extracted from the 
trajectories according to /S.E = FdragdL' where Lq is 
the reference point and L is an arbitrary point along the 
trajectory. We used the Stokes’ drag formula for a sphere 
equally immersed in the subphase fluids, F^rag = bTr/xUa, 
where /i is the average viscosity of oil and water evaluated 
at the temperature of the environment. The curvature 
capillary energy found over the trajectories beginning at 
aAco = 6 X 10“^ and ending ten particle radii from the 
micropost was plotted against aAco as open symbols in 
Fig. [3 this range of aAco was selected because all trajec¬ 
tories captured in experiment spanned this region. The 
relationship is linear and the total energy difference along 
a typical trajectory is thousands of times greater than 
thermal unit energy ksT. The curvature capillary energy 
is presented normalized by qTra^ = 8.8 x 10^ ksT] the en¬ 
ergy for the segment of the particle paths shown in Fig.[^ 
is in the range of 6, 000 — 50, 000 ksT. This magnitude in¬ 
dicates that the equilibrium wetting boundary condition 
cannot be responsible for the migration of the spheres 
in our experiment, as according to Eq. for a typical 
microsphere a = 5/xm with equilibrium contact angle of 
0 = 90°, the curvature capillary energy is of magnitude 
/S.F ^ sin^ Oq A4 lO/c^T. Moreover the relation¬ 
ship between /S.E and roAco would be highly nonlinear. 
Hence, we conclude that the equilibrium contact angle 
boundary condition does not apply to our microspheres. 
Rather, the microspheres migrate with pinned contact 
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FIG. 8. (a) Time-stamped trajectories of microspheres 

around a square post illustrating that the microspheres fol¬ 
low complex trajectories as defined by the deviatoric curva¬ 
ture field (the scale bar is lOO/xm). Numerically evaluated 
curvature gradient around (b) a square post and (c) a corner 
of the square post. The arrows scaled with the magnitude of 
deviatoric curvature gradient. Movies of particle migration 
are provided as supplementary material. 


lines. 

We propose that this curvature migration is an as¬ 
say for nanometric corrugations of the contact line in a 
trapped state. The magnitude of the quadrupolar mode 
hp for the contact line undulations can be inferred from 
the trajectories in Fig. While 7 of 10 of the reported 
trajectories have hp similar in scale to the particle rough¬ 
ness (between 20 — 40 nm), magnitudes for the remaining 
trajectories were larger, with hp as high as 130 nm for 
one trajectory. These results indicate that similar par¬ 
ticles have differing pinning states at the interface, with 
significant consequences for their ensuing dynamics. The 
role of gravity for particles on curved interfaces has been 
addressed previously [HUS- Because the Bond number 
is negligible, particle weight plays no role in the defor¬ 
mation of the interface, and the analyses above are valid. 
For weak enough curvature gradients, however, particles 
cannot overcome potential energy barriers and thus can 
attain an equilibrium height. 

This form for the curvature capillary energy has been 


invoked previously to study curvature capillary migra¬ 
tion of cylindrical microparticles which followed complex 
trajectories on interfaces around square microposts with 
associated complex curvature fields [2]. If the spheres 
indeed have identical physics, they, too, should migrate 
along complex contours in such a curvature field. We 
have studied trajectories of spheres in this setting (see 
Fig. [^a), in which isolated particles migrate to corners, 
as does a pair of dimerized particles at the lower right 
hand corner). To compare particle trajectories to lo¬ 
cal curvatures, the curvature field around the square mi¬ 
cropost was determined using a Galerkin finite element 
method, as discussed in Appendix C. Vectors indicating 
the magnitude and direction of gradients in deviatoric 
curvature are indicated in Fig.j^b) and (c). The spheres 
migrate from their initial point of attachment along paths 
defined indeed by these vectors. The ability of the sphere 
to trace this complex trajectory confirms the underly¬ 
ing physics of microparticles with pinned contact lines is 
similar regardless of details of the particle shape, and is 
consistent with the concept that the particle quadrupolar 
mode couples to the underlying saddle shaped surface. 


V. CONCLUSION 

We study micropshere migration owing to curvature 
capillary energy in theory and experiment. We show that 
for equilibrium contact angles, the curvature capillary 
energy is very weak, with leading order contributions of 
fourth order in deviatoric curvature or higher, in contra¬ 
diction to the accepted form in the literature. This lead¬ 
ing order contribution would amount to roughly ksT 
in our experiments, several orders of magnitude lower 
that we measure in experiment. In experiment, micro¬ 
spheres migrate along deterministic trajectories defined 
by curvature gradients. We find that the corresponding 
capillary curvature energy propelling the particles ranges 
from 6,000 — 50,000 ksT. We compare these observa¬ 
tions to arguments derived previously for particles with 
pinned contact lines, in which the quadrupolar mode of 
the particle contact line undulation couples with the cur¬ 
vature field to yield an energy linear in the deviatoric 
curvature. The data indeed obey this form, allowing the 
magnitude of the particle induced quadrupole to be in¬ 
ferred. In many cases, it is comparable to particle rough¬ 
ness as determined by AFM. However, significantly larger 
magnitudes are also found, suggesting that similar par¬ 
ticles can have different pinned states at the interface. 
These results imply that contact line pinning occurs for 
microparticles at these curved fluid interfaces with dra¬ 
matic implications in their dynamics at interfaces. 


VI. APPENDIX A 

In this section, we provide a detailed derivation for 
the interface shape when a particle of radius a which at- 














FIG. 9. Schematic representation of inner y, z) and outer 
(X, V, Z) coordinates. 


tains an equilibrium contact angle is placed on a host 
interface. The magnitude of this disturbance created 
by the particle is comparable to the particle radius, a. 
The disturbance is a monotonically decaying function 
which decays over distances comparable to a. More¬ 
over, the particles are small compared to the character¬ 
istic length of the host interface Rc ^ Vacq and hence 
the interface domain can be divided into two distinct re¬ 
gions, an inner region close to the particle where the rel¬ 
evant length scale is a and an outer region further away 
where the relevant length scale is Rc. For small ratio of 
A = = roAco ^ 1, the problem can be solved sys¬ 

tematically by a matched asymptotic expansion scheme 
[Hilo]. We perform this matching here, in detail, for 
our particular interface shape: In the outer region (far 
from the particle) the shape of interface is known and 
the shape is given according to a logarithmic expression, 

j^outer ^ ^ ln(, (33) 

Rc Rc Rm 

please note that the above relation is given in the (outer) 
coordinate with origin at the center of the circular mi¬ 
cropost (X,Y, Z) where Z = and they are nondi- 

mensionalized with respect to the host interface radius of 
curvature Rc. 

In the inner region on the other hand, we adopt inner 
coordinate (x, y, z) defined with respect to the tangent 
plane which has slope of e = _Rrn^^}X ^i^h respect to 
the outer coordinate (See Fig. and Lq is the location 
of particle center of mass with respect to the center of 
micropost.. This coordinates is made dimensionless with 
inner region length scale ro = asin^o- The height of 
the interface in the presence of the particle is defined 
~ _ j^inner^ where we utilize a Monge representation of 
the interface. The shape of interface in the inner region 
can be governed via Young-Laplace equation which, in 
the limit of small gradients and infinitesimal Bond num¬ 
bers, reduces to: 


(34) 


We then expand the non-dimensional height of the inter¬ 
face in the inner region as a power series to yield, 

= ho + hiX + h2X^+ 0{X^). (35) 

The equilibrium wetting boundary condition at three 
phase contact line is dictated by Young’s equation, 

np ■ n/= cos^ol^gp, (36) 

where np is the unit normal to the particle and n/ is the 
unit normal to the interface which can be determined as. 


np = ep = sin Obc + cos Obz 


n/ = 


^ s h 


(37) 

(38) 


where (i?, 0, 0) and (r, 0, z) are spherical and cylindrical 
coordinates located at the particle center of mass, respec¬ 
tively. The operator is the surface gradient defined 
as. 


V, = (/-e,e,).V, (39) 

and / is the unit tensor. By replacing Eq. and in 
Eq. |37| and |38l Eq. |36| reduces to. 


I _cos0o + ^-fe,-V,/i(r)| 

np • n/Irgp -- / -lrep> 




cos00 + sin0o[^ - ■ Vs/i(r)] 


, I, dh{v 


r dcj) 


+ (40) 


fe(r) gfe(r) 


ro 


ro dr 


where we have utilized the small gradient approximation 
|Vs/i(r G p)| <C 1. Einally, the boundary condition (in 
dimensional form) at particle contact line can be written 
as. 


cos^0r/5^x9 


sin 0o\ — 
ro 


(1 — cot 0 


hdh 

ro or 


(41) 


By recasting this expression in dimensionless form and 
substituting it in Eq. we can collect terms of similar 
power in A. One can show. 


h 

h 


dhi 

dr 

dh2 

dr 


r=l 


f=l 

cot Of) ,,dhi^ A dhi 


2 

) -2hi 


dhi 

dr 


(42) 


y2'^inner ^ 


f=l 


(43) 
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where Eq. 42 was derived by Wiirger m, and is given 
in Section 2 in dimensional form as Eq. the bound¬ 
ary condition at the three phase contact line. Eq. is 
the first correction to this contact line boundary condi¬ 
tion. In the above, we assumed Jiq = 0. This hypoth¬ 
esis is proven below. The far field boundary condition 
for the inner solution is obtained by matching the inner 
solution with the outer solution via Van Dyke’s match¬ 
ing procedure [191 ED]- This can be done systematically 
by writing the outer solution in terms of inner variable; 
thereafter holding the inner variable fixed and expand¬ 
ing for small parameter A. This provides the far field 
boundary condition for the inner region. Eormally, this 
matching is written: 


lim = lim h 

r^oo A^O 


outer 


(44) 


note that in the above expression, is made dimen- 

sionless with outer region length scales, and is 

made dimensionless with inner region length scales. 

We demonstrate the execution of this matching proce¬ 
dure for our particular host interface, without loss of 
generality. The same procedure can be applied to any 
host interface provided ^ ^ 1. To perform this match¬ 
ing, we perform a coordinate transformation between the 
inner and the outer coordinates. The inner coordinate is 
tangent with respect to the outer coordinate with slope e. 
The relationship between the two coordinates providing 
small slopes |e| ^ 1 (in dimensional form) can be written 
as. 


X — Lq X 0(e), 

Y = y, (45) 

Z = Zq z 0(e), 

where Zq is the interface height at the particle center of 
mass with respect to the outer coordinate. Utilizing the 
above transformation, the right hand side of matching 
condition expression in Eq. [^can be written as. 


implies exploring regions around the particle large com¬ 
pared to the particle radius. To apply the boundary con¬ 
dition, we introduce a vertical shift factor ujq = ^ 
to the particle center of mass and therefore we have. 




cos 20 
4 


(f 2 + 



cot 00 cos 4^ 

~36 


(49) 

(50) 


This is the inner solution to the shape of the interface up 
to A^. 

To construct a uniformly valid (uv) solution over the en¬ 
tire domain, we define a uniformly valid solution in di¬ 
mensional form according to, 

^ - Rc lim , (51) 

A^O 
f fixed 

consequently, the (dimensional) disturbance will be, 


n = h'^'' - R^h 


outer 


= A 


^0 2^ocot6>o 


12 f 2 


cos 20 — A 


36f4 


cos 40, 
(52) 


the disturbance is solely a decaying function of inner vari¬ 
able f. This implies that its value is identically zero in 
the outer region. Thus, the particle results in a “local” 
disturbance which fades over a length scale comparable 
to the inner length scale tq. Eurthermore, in evaluating 
the area owing to the particle disturbance Eq. there 
is no coupling between and owing to orthogonal¬ 
ity of cosn0. Therefore, there is no term in the energy 
expression of order A^. In the text, we evaluate the en¬ 
ergy up to A^ and find this term to have prefactor zero. 
Hence, the leading order contribution to the curvature 
capillary energy for spherical particles with equilibrium 
contact lines is order A^. 


VII. APPENDIX B 


Re 


i^^tan^ ,v(^ + To) 
■ ln( — 


Re 


Ru 


Re' 

(46) 


By nondimensionalizing and expanding the argument of 
the log term for small A, one can show. 


According to the asymptotic discussion, all area in¬ 
tegrals for the free energy can be decomposed to two 
contributions in the inner and outer regions: 



/ 


inner 


dA + 


fouterdA_ ( 53 ) 


inner 


outer 


lim 0) = cos20 + 0(e,A^), (47) 

A^o 4 

r fixed 

utilizing Eq. this can be recast to specify the far field 
boundary condition for the disturbance in the inner re¬ 
gion, 

lim (r, (p) = cos 20 + 0(e,A^), (48) 

r^oo 4 

in dimensional form, this becomes Eq.[^ This matching 
also clarifies the meaning of taking the limit r ^ oo. This 


In particular the integral which accounts for the inter¬ 
action of the disturbance with the host interface can be 
written: 

Vho • XrjdA = 

D-P 

(V/io • V7?)*"""’'dA -h ^ (V/io • VriY'^^^^dA, (54) 

since, as shown in Appendix A, the disturbance de¬ 
cays to zero in the inner region, the second integral 
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in the right hand side of Eq. is zero. Asymptotics 
also establish that the description of the host interface, 
lim ho = ^ cos 20 remains valid ev- 

erywhere that the disturbance is finite. This allows the 
first integral on the right hand side of Eq. be evalu¬ 
ated. In the ensuing discussion, the superscript “inner” 
is removed, while understanding that the sole surviving 
integral is in fact in the inner region. Substituting ho 
and 7^, the integral can be evaluated directly. Or, the 
divergence theorem can be used to recast this integral in 
terms of contour integrals that bound the domain; there 
are two equivalent forms for this approach. We evaluate 
the integral all three ways: 

(i) Direct integration over the angle 0 shows this inte¬ 
gral to be zero at all radial locations between ro and any 
location r* in the domain. 


which, via application of Green’s theorem as: 
V • {hoVr])dA 


D-P 


= lim (p Br • (hoV7^)rd0 — 0 • {hoVr])rd(l)^ 

r*^oo J J 


(59) 


the two contour integrals are: 


lim / Br ' {r]S/hQ)rd(j) = — , (60) 

J 24 


-/ 


• hoWrjrdcj) = 


24 

7rro"^Aco^ 

24 


, (61) 


nr* n2'K 

Jro Jo 


dr dr ^ r2 dd) ^ 


rr* A 2 4 p27T 

/ —^——dr / (—cos^20 + sin^20)d0 

Jro 12 r Jo 


rn'^Acn^ r* 

^ln(-)(0)=0. 


(55) 


Thus, for any slice of the interface around the particle, 
this integral is zero. One might at first be concerned 
about the logarithmic divergence in Eq.[^ but the factor 
multiplying it is identically zero for any value of r*. 

(ii) One form of integration by the Green’s theorem. The 
integrand can be recast as: Vho • = V • {rj'Vho) — 

( 77 ^ 7 ^ho), where the second term on the right hand side 
is zero since V^ho = 0. The area integral can then be 
recast: 


Vho • VrjdA = ® V • {r]Vho)dA, (56) 


D-P 


D-P 


the integrals can be written via application of Green’s 
theorem as contour integrals: 


V • {r]\/ho)dA 

D-P 

= lim 0 • {r]Vho)r d(j) — 0 • {r]Vho)r d0, 

r*^oo J J 

r=r* r=ro 

(57) 


which ultimately yields the results given in Eq. and 

ES 

(iii) Another form of integration by Green’s theorem. 
The integrand can also be recast as Vho • V 77 = V • 
(hoVr/) — (hoV^? 7 ), where the second term on the right 
hand sides is zero since V^rj = 0. The area integral is 
then. 


Vho • VrjdA = 


V • {ho\/r])dA, (58) 


The terms are equal and opposite and hence the integral 
is zero. In the literature, this area integral was evaluated 


by first recasting it in contour integrals as in Eq. 57 


Thereafter, the far field contour was neglected, and the 
quadratic term in Acq was retained as the value for the 
integral (see Eq. 23). If the same logic is applied to the 
area integral as recast in Eq. only Equation Eq. 
would be retained. These two forms, which should be 
equivalent, are in fact equal in magnitude and opposite 
in sign. This clearly demonstrates that the neglect of the 
outer contour is improper, as it represents a failure to 
close the bounding contours on the integral, and violates 
the Green’s theorem used to recast the area integral. 


VIII. APPENDIX C 


In this section we discuss the numerical method used to 
calculate the height of interface around the square post. 
The height of interface in this case can be governed via 
the linearized Young-Laplace equation according to, 

= Apgh + Ap, (62) 

where Ap is the density difference of subphase fluids and 
Ap is the pressure jump across the interface. We scaled 
the variables as follows, 

" = 1M- ^ = 

Apg 

Eq. I^can be rewritten as, 

V^h-h = ±l, (64) 

where ±1 on the right hand side stands for negative and 
positive pressure jumps across the interface. To obtain 
the discretized form of Eq. |^in the weak form to imple¬ 
ment the finite element algorithm for solution, we first 
multiply both sides of the Eq.j^by a weighting function 
(fi and then integrate over the area of element. 


(piV‘^hdA ■ 


(fihdA = ± 



(65) 


D-P 


D-P 












11 



FIG. 10. Numerical solution of the curvature field via finite 
element analysis around a square post, (a) Deviatoric curva¬ 
ture and (b) mean curvature. The magnitudes of the devia¬ 
toric and mean curvatures are in . 


Having utilized integration by parts and the divergence 
theorem, one can show, 


(piVh • mds — (h) Vifi • VhdA — (h) ipihdA = 





( 66 ) 


where m is the unit normal vector of the enclosing con¬ 
tours which are pointing outward from the domain of in¬ 
terest. In the above expression, the first term only applies 
to the boundaries and therefore we neglect it since we 
have Dirichlet boundary conditions in our problem. Uti¬ 
lizing the Galerkin method [25j , the non-dimensionalized 
height h were approximated by a piecewise quadratic 
polynomial function where = 1 at vertex i and ipi ={) 

at all other vertices. Thus, the non-dimensional height 

N 

can be approximated as where N is the 

1=1 

total number of vertices in the mesh and hj is the value 
of the height at the jth vertex. Since the problem is 2D, 
we generated a triangular structured mesh by discretiz¬ 
ing the domain. This was done in such a way that more 
nodes were available in the transition region (sharp cor¬ 
ner) to capture the rapid change in height. The pressure 


jump was not known a priori; it is found from the slope 
condition at the middle of the square post. To do so, we 
solve the entire problem in an iterative manner. We guess 
a value for the pressure jump across the interface, solve 
the problem by applying the Dirichlet boundary condi¬ 
tions at the post and the outer ring, and subsequently 
calculate the slope at the middle of the square post. The 
slope is compared to the measured value in the experi¬ 
ment with ip = 15°. The new pressure jump is estimated 
according to the following, 

(67) 

midplane 

We follow this procedure until we reach the convergence 
and the difference between two consecutive iterations be¬ 
comes less than 1%. The pressure jump we found is 
Ap = —0.2170^. The solution was also tested against 
the number of triangular meshes to confirm that the final 
result is independent of number of meshes in the domain 
of interest. Furthermore, to evaluate the principle cur¬ 
vatures, we evaluated higher order partial derivatives in 
each individual mesh according to, 

^ o2 
hxx = ^ 
i=l 

Thereafter, we calculate the principle curvature accord¬ 
ing to appropriate relations. The resulting deviatoric and 
mean curvature fields are shown in Fig. the devia¬ 
toric curvature diverges near the corners, spanning val¬ 
ues ranging from 400 m~^ to 3000 m~^ whereas mean 
curvature varies from 3 m~^ to 5 m~^. It is therefore 
safe to ignore capillary energies owing to mean curvature 
gradients. 
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